NASA Technical Memorandum 101482 


Finite Element Applications to Explore the 
Effects of Partial Bonding on Metal 
Matrix Composite Properties 

< l&Zl - It - 1C 14 t 2) f i O'] £ tic till hb'i-i.Ctlt 

HtUtAlKli r iC tJiiCfi Uf HitC'M OF 

l Ibn PL I C Mi ft C t filH P* it i) cC f tCSlli 

UCItflJM (MiA) 26 p COLL 1 1 L* Gr.cld£ 

G3/24 ClbSCeC 


J.J. Caruso 

Lewis Research Center 
Cleveland, Ohio 


D. Trowbridge 
University of Akron 
Akron, Ohio 


and 


C.C. Chamis 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

30th Structures, Structural Dynamics and Materials Conference 
cosponsored by the AIAA, ASME, ASCE, AHS, and ACS 
Mobile, Alabama, April 3-5, 1989 




E-4610 


FINITE ELEMENT APPLICATIONS TO EXPLORE THE EFFECTS OF PARTIAL BONDING ON 
METAL MATRIX COMPOSITE PROPERTIES 

3 . 3 . Caruso 

National Aeronautics and Space Administration 
Lewi s Research Center 
Cleveland, Ohio 44135 

D. Trowbridge 
The University of Akron 
Akron, Ohio 44325 

and 

C.C. Chamis 

National Aeronautics and Space Administration 
Lewi s Research Center 
Cleveland, Ohio 44135 


Abstract 

The mechanics of materials approach (defini- 
tion of E, G, v, and a) and the finite element 
method are used to explore the effects of partial 
bonding and fiber fracture on the behavior of high 
temperature metal matrix composites. Composite 
ply properties are calculated for various degrees 
of disbonding to evaluate the sensitivity of these 
properties to the presence of fiber/matrix disbond- 
ing and fiber fracture. The mechanics of mate- 
rials approach allows for the determination of the 
basic ply material properties needed for design/ 
analysis of composites. The finite element method 
provides the necessary structural response (forces 
and displacements) for the mechanics of materials 
equations. Results show that disbonding of frac- 
tured fibers affect only Elll and aft.] ] 
significantly. 

Nomenclature 

A area 

E Young' s Modulus 

e effective value 

F force 

f*m,Q, fiber, matrix, ply 
G shear modulus 

S physical dimension of the model 

t temperature 

u,v,w displacement in 1, 2, 3, direction 
X,Y,Z global coordinate axis system 

a thermal expans ion . coeffi ci ent 

c strain 

v Poisson’s ratio 

a stress 

1,2,3 Laminate coordinate axis system 


with finite element analysis to determine the 
behavior of epoxy matrix composites J This method 
has been extended to high temperature metal matrix 
composites. 2 

The finite element method is not intended to 
replace experiments; but* to provide needed 
information not yet available through experiments. 
The finite element method provides a means to 
computationally simulate a physical experiment to 
estimate quantities which are difficult, or 
impossible, to menasure. 

The objective of this work is to study the 
influence of disbonding of fractured fibers on the 
ply properties of metal matrix composites. In 
the study the disbonded fiber is considered frac- 
tured and therefore the load transfer between the 
fiber and matrix is inhibited by the amount of 
disbonding. 

Model i ng 

The composite system considered for this work 
is PlOO-Graphi te/Copper at 0.466 fiber volume 
ratio. Analyses for high temperature and room 
temperature ply material properties are conducted 
using Version 65 of MSC/NASTRAN . 

A useful analytical tool in the evaluation of 
structures with repeated geometry is finite element 
analysis using the superelement method. 2 This 
involves the partitioning of a finite element mesh 
into separate collections of elements called 
superelements. Each superelement is solved sep- 
arately and then combined to complete the 
analysis . 

The composite superelement mesh is of modular 
nature, consisting of a unit cell (Fig. 1) and 
images of the unit cell; therefore, discontinu- 
ities can easily be added by substituting a 
conventional mesh, which models the di sconti nui ty , 
in place of one of the images. This conventional 
mesh can also be imaged to simulate multiple 
discontinuities within the composite material. 

Some discontinuities of interest that can be 
modeled in this manner are partial bonding and 
fiber fracture, both common in metal matrix 
composite materials. 

The cases chosen are but a few of a multitude 
of possible partial bonding configurations. 
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However, the cases selected are sufficient to 
understand the effects of partial bonding on the 
composite's material properties. 

The mesh that is employed in this work con- 
sisted of a cluster of nine fibers in a three by 
three matrix (Fig. 2). The initial version of 
the mesh is assembled with the center fiber/ 
matrix block modeled as a primary superelement 
and all surrounding blocks are images of the 
primary. This mesh is used when no fibers are 
disbonded and when all nine fibers are simulta- 
neously disbonded. The mesh is reassembled when 
only the center fiber of the nine cell model is 
disbonded. For this case, the center block is 
modeled using conventional finite elements. The 
primary superelement is defined to be to one side 
of the center fiber and the remaining perimeter 
superelements are modeled as images of this one. 
Minor modifications are made to each of the meshes 
to allow for the modeling of varied amounts of 
disbond between fiber and matrix. 

The approach used to simulate the disbond is 
described briefly here. Duplicate grid points are 
defined for all grids around the circumference of 
the fiber, at the interface between fiber and 
matrix. This duplicate grid point is included in 
the analysis only when the portion of fiber where 
it resided is to be disbonded from the matrix. 

The disbonding took place when the connectivity 
cards that define the perimeter of the fiber are 
altered so that one of the duplicate grids is 
associated with the fiber, and one with the 
matrix. No connectivity then exists across the 
interface of constituents, in effect producing a 
crack of zero width, with fiber on one side and 
matrix on the other. In the work described here, 
the total ci rcumference of a fiber is released 
together. The amount of disbond is varied by 
disbonding different lengths of fiber as described 
above. For example, 2.78 percent disbonding of 
the center fiber is detailed in Fig. 3 with the 
plane of symmetry where the grids are not dis- 
bonded. No attempt is made to find the loading at 
which the disbonding or fiber fracture will take 
place or to prevent the overlapping of the fiber 
and the matrix after disbonding occurred. Instead, 
interest is placed on the effect of disbonding of 
the fractured fiber on the composite's material 
properties . 

Procedure 


INTRODUCTION 

Recent research at NASA Lewis has led to the 
development of a computational method to predict 
the behavior of unidirectional composites. The 
method utilizes a mechanics of materials approach 


where si is the length in the 11 direction. 


2 is found by finding the average deflec- 
tion in the Y direction as a result of the 
enforced displacement in the X. Dividing this 
deflection by the width of the model yields a 
strain in the transverse direction (22). The 
Poisson's ratio v<i\2 15 calculated from: 


v ai3 


V 


ai2 


e Q.22 
e 41 1 


is found by similar methods. 


( 2 ) 


To calculate the Poisson's ratio and 

the transverse modulus Eg, 22 * the f ace w i th the 
negative Y direction as its normal is fixed in 
the Y direction (v a = 0.0) and the face with the 
positive Y direction as its normal has an 
enforced displacement in the positive Y direction 
<Vg = v) resulting in tension transverse to the 
fiber direction (Fig. 2). The equivalent applied 
stress is then calculated by averaging the result- 
ant forces over the face and dividing by the 
face's area. E ^22 and are calculated 

similarly to Eqj ] and vju 2 respectively. 

E &33 and vq ( 32 are calculated from another 
unique set of boundary conditions in the same way 
E 0,2 2 and v £23 are calculated. 

For the determination of G&21 » an enforced 
displacement in the X direction is placed on the 
face with the positive Y axis as its normal. The 
shear strain is then calculated by dividing this 
displacement by the width of the model. The 
effective shear stress in the 21 direction is then 
calculated by dividing this force by the area over 
which it is applied. The shear strain for small 
displacements is given by the enforced deflection 
divided by the width of the specimen. G^ 2 l is 
obtained by dividing the shear stress in the 21 
direction by the shear strain. 


The following are the procedures used to cal- 
culate ply material properties from the finite 
element output. The subscripts "11" are used to 
denote the direction along the fiber (longitudinal) 
and the subscripts "22" and "33" are used to 
denote the directions transverse to the fiber. 

vqj 3 , VQJ 2 > and Eqji are obtained from the 
same loading conditions. In this case the front 
face is fixed in the X(u] = 0 . 0 ) and the back 
face is displaced in the X(ug * u) (Fig. 2). 

Eqji is calculated by finding the total 
force over the displaced face and dividing this 
force by the area of this face yielding an equiv- 
alent applied stress. The strain is then cal- 
culated by dividing the applied displacement by 
the length of the specimen. Young's modulus is 
computed from the mechanics of materials equation 
in the longitudinal 11 direction by: 


G 121 


e 

°U\ 

c ii21 


G& 3 i is found by similar methods from 
another unique set of boundary conditions. 


(3) 


G£23 can be obtained by applying an enforced 
displacement in the positive Z direction to the 
side whose normal is the positive Y axis. The 
total resulting force on this face is calculated 
from the finite element output. This total shear 
force in the 23 direction is divided by the area 
of the side over which it is applied, resulting 
in the average shear stress in the 23 direction. 
The small displacement shear strain in the 23 is 
then calculated by dividing the applied displace- 
ment in the Z direction by the width of the 
model. Gq ,23 is calculated by dividing the 

effective shear stress, a^ 23 , by 
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In order to calculate the thermal expansion 
coefficients a^n , ag, 22 . dn( * a a 33 > a P^ ane 
the center of the model, with the positive X 
direction acting as its normal, is restricted 
from movement in the X direction (U5 = 0.0). 
This restriction assures symmetry about the 
center of the model. A thermal load is then 
applied to the model, T = T 0 . 

The first step in finding agjj is deter- 
mining the average of the displacements on the 
end with the positive X axis as its normal and 
dividing by the length of the model, S], 
resulting in the strain in the 11 due to the 
change in temperature. Now, by dividing the 
strain by the change in temperature, one can 
obtain the thermal expansion coefficient in the 
11 direction. 


*122 is found by taking the average dis- 
placement in the Y direction and dividing by the 
change in temperature and the width. 

The calculation of 0^33 is similar to that 
of a£ 22 • with the exception being that the aver- 
age displacement in the Z is divided by the 
change in temperature and by the height of the 
model . 

A set of reference finite element runs for 
PlOO/copper are executed with a 0.466 fiber volume 
ratio. These runs contained no disbonding or 
fiber fracture. Composite ply material properties 
are calculated from the output of these finite 
element runs (Table I). Additional confidence was 
gained by comparing the ply properties in Table I 
with those predicted by Caruso and Chamis 2 and by 
Hopkins and Chamis 4 with very good correlation. 

Disbonding and fiber fracture of a single 
fiber in the center of the nine cell model is now 
considered (Fig 3). Room temperature (70 °F) 
material properties are used (Table II). In all 
loading conditions, fibers that are to be dis- 
bonded are also considered fractured and; there- 
fore, did not have an applied displacement 
assigned to them. This modeled the fiber fracture 
in advance of disbonding. In this first series 
of runs, each element around the circumference of 
the fiber is allowed to disbond using the methods 
described earlier. The disbonding of a fractured 
fiber is done for each successive layer of ele- 
ments along the length of the fiber. After the 
mesh is changed to reflect the next layer of dis- 
bonding, the loading conditions are applied and 
ply material properties calculated. Results are 
plotted for each of the ply material properties 
verses percent disbond of the total mesh (Figs. 4 
to 7). These plots show the properties affected 
most by disbonding and to what degree they are 
affected. Some of these properties can be used 
as indicators of partial bonding of a fractured 
fiber in cases where material properties are not 
near predicted values. When suspecting the 
presence of partial bonding, closer attention 
should be paid to those properties found to be 
more sensitive to the presence of the fiber 
disbonding. 


The next series of runs are executed with 
high temperature constituent material properties 
(Table II). The temperature selected is 1500 °F. 

The results of the material property calculations 
are shown in Figs. 8 to 1 1 . Note that these plots 
reflect a little over 11 percent disbonding of the 
entire composite. This is a result of releasing 
only the center fiber for these cases considered 
here. The 11 percent disbonding value shown on 
these plots represents the point at which the 
center fiber is totally disbonded and the surround- 
ing eight fibers are fully bonded to the matrix. 

When the fractured center fiber is considered 
completely disbonded it is still connected to the 
matrix material by a ring of nodes in the plane 
of symmetry (Fig 3). This can add some stiffness 
transversely and in shear. 

The next series of conditions considered are 
when all nine fibers in the model are disbonded 
and the room temperature boundary conditions are 
imposed. In this case the percent disbonding can 
go up to 100 percent of the total fiber length. 

Once again, the connectivity of a ring of nodes in 
the plane of symmetry around each fiber is main- 
tained at 100 percent disbonding. Material pro- 
perties are calculated by applying each of the 
loading conditions described earlier. These 
results are plotted showing the effect of fiber 
disbonding on each property when all fibers are 
equally disbonded (Figs. 12 to 15). 

The next series of runs is repeated for high 
temperature conditions. All nine fibers are 
released as described above using material proper- 
ties that reflected a use temperature of 1500 °F. 
Material properties are calculated and plotted for 
varying degrees of fiber disbonding ranging from 
0 to 100 percent disbonding (Figs. 16 to 19). 

Results and Pi scussion 

The degradation in the ply material proper- 
ties could be influenced by the difference in the 
corresponding fiber and matrix material proper- 
ties. As the fibers become structurally inactive, 
the matrix material properties have a more signi- 
ficant role in determining the composite's global 
properties. Ideally, the composite's material 
properties are bounded by the highest and lowest 
values of its constituents; however, it will be 
seen that a composite ply material property can 
be lower than the lowest constituent property. 

Room Temperature Ply Properties with Center Fiber 
Pi sbondinq 

Figure 4 shows the degradation of the modulus 
in the 11 , 22 and 33 directions due to disbonding 
of the center fiber only. These results reveal 
that the longitudinal modulus has decreased by 
about 8 percent while disbonding is at 11.1 percent 
of the total circumferential fiber area of the 
model . 

The percent decrease in E&22 and E £33 
for 11.1 percent total fiber area disbonding is 
1.5 percent. From these results it can be con- 
cluded that E an is more sensitive to fiber dis- 
bonding than the transverse moduli. The decrease 
in E^n is due mostly to geometric effect. The 
geometric effect is the holes left behind by the 
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disbonded fiber much like a sponge. From further 
results it will become evident that. E^n is the 
second most sensitive property to fiber 
disbonding. 

The effect of fiber disbonding on Poisson's 
ratio for center fiber disbonding is depicted in 
Fig. 5. When the center fiber is completely 
disbonded ( 11.1 percent total fiber length dis- 
bonding), v &]2 and vq > ]3 have degraded by 
0.5 percent. A decrease of 1.6 percent in vq. 23 
and v £32 is witnessed when total center fiber 
disbonding is achieved. This degradation of 
Poisson's ratio is very small and can be con- 
sidered insignificant. 

The decrease in shear modulus due to center 
fiber disbonding is shown in Fig. 6 . The percent 
change in G &12 is 3 percent and G &23 is 
2 percent when the center fiber is completely 
disbonded. This shows that G&| 2 , 0^3 as well 
as g &23 are relatively insensitive to the dis- 
bonding of the center fractured fiber. 

The thermal expansion coefficients are con- 
sidered in Fig. 7. The lines representing the 
transverse thermal expansion coefficients are 
coincident, as to be expected. The disbonding 
conditions considered here produced a decrease in 
the transverse thermal expansion coefficients by 
0.3 percent. This decrease is minimal and may be 
considered to be insignificant. 

The line describing the change in the longi- 
tudinal thermal expansion coefficient is also 
displayed in Fig. 7. ccqj i increased by 20 per- 
cent when the center fiber is totally disbonded. 

Of all the material properties considered, a&n 
is the most sensitive to center fiber disbonding 
and i s very significant. 

High Temperature (1500 °F) Ply Properties with 
Center Fi ber Pi sbondi nq 

Figure 8 shows the degradation of longitu- 
dinal and transverse moduli. Since 11.1 percent 
of the effect of the fiber has been removed due 
to the total disbonding of the center fiber, a 
decrease of only a 7.3 percent is witnessed in 
E <i}}. This decrease is of the same magnitude as 
that observed for the room temperature E^i i . 

The conclusion is that, the change in Eq.ii 
brought about by center fiber disbonding is 
significant. 

Due to geometric effects, Eq .22 decreases by 
2 percent when the center fiber is disbonded 
(Fig. 8 ). This decrease is due to the reduced 
effective cross-sectional area, and the subsequent 
reduced transverse stiffness of the composite, as 
resulting from the disbonding of the center fiber. 
This hole removes some of the stiffness from the 
model resulting in a lower transverse modulus. 

Note that this decrease is almost the same as the 
room temperature decrease in Eq .22 and £ 0 , 33 - 

In Fig. 9 the Poisson's ratios of the com- 
posite are plotted, vqj 2 and vq .23 decrease by 
approximately 1.5 percent when the center fiber 
is disbonded fully. When compared to the room 
temperature results, it is evident that there is 
a small increase in the change of V£| 2 > and the 
change in vq. 23 is about the same. From these 


results it may be concluded that temperature has 
little effect on the rate of degradation of 
Poisson's ratios, vg, ] 2 » v &13* v &23 and ^ 0,32 
may all be considered to be relatively insen- 
sitive to center fiber disbonding. 

G &| 2 degrades by almost 5 percent with the 
disbonding of the center fiber (Fig. 10). The 
change in G^l 2 at r00m temperature is 3 percent 
and at high temperature it is 5 percent when the 
center fiber is disbonded. G^2 3 * at h ^ 9 h temper- 
ature, has degraded by almost 3 percent when the 
center fiber has totally disbonded. Gqj 2 an d 
Gq, 23 can de said to be relatively insensitive 
to center fiber disbonding. 

Figure 11 contains the plots that represent 
the changes in thermal expansion coefficients at 
high temperature as the center fiber is disbonded. 
a £22 increased by less than 1 percent, an insigni- 
f i cant amount. 

In figure 11, ot&i] shows an increase of 
17.3 percent when the center fiber is completely 
disbonded which represents a decrease in effective 
fiber area of 11.1 percent. This change in 
a^ii is significant and close to the percent 
change in ctQj 1 at room temperature. 

Room Temperature Ply Properties with All Fibers 
Pi sbondi nq 

The effect on the ply moduli of all the 
fibers disbonding is plotted in Fig. 12. The 
degradation of Eq,h is now obviously nonlinear. 

Efc] 1 , with 100 percent fiber disbonding, has 
decreased by 78.2 percent and is lower than E m . 

The cylindrical holes left by the disbonded fibers 
decrease the stiffness of the composite, much like 
a sponge. The composite can actually have a lower 
modulus than that of the matrix alone. 

E ^22 decreased by only 15 percent with 
100 percent fiber disbonding (Fig. 12). The 
decrease in E ^22 is due *° the sponge effect. 

In Fig. 13 vgj2 increases until about 
25 percent disbonding, then VQ .12 decreases grad- 
ually to 100 percent disbonding. This phenomenon 
will have to be examined in future work. 

One hundred percent disbonding brought about 
a decrease in v&23 of 26.7 percent as shown in 
Fig. 13. The value at 100 percent disbonding is 
less than that of the matrix alone. This may be 
due to the spongy effect of the disbonded com- 
posite, where a significant percentage of the 
transverse displacements take place in the cylin- 
ders left by the disbonded fibers, resulting in 
less transverse strain. 

Figure 14 represents the results of disbond- 
ing on the shear moduli at room temperature. Both 
G &}2 and 6^3 show an approximate degradation of 
25 percent when the fibers are completely disbonded. 
As before, at 100 percent disbonding, these ply 
properties are lower in value than that of the 
matrix material alone, presumably due to the 
sponge effect discussed earlier. 

Figure 15 shows the change in the composite's 
thermal expansion coefficients when disbonding is 
taken from 0 to 100 percent at room temperature, 
aim increases drastically as the disbonding 
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increases. At 100 percent disbonding, a&]] has 
increased by 557 percent. The value of a^}] at 
100 percent disbonding is virtually identical to 
that of the matrix alone. 

The degradation of <x&22 a ^ so depicted 
in Fig. 15. a&22 degrades by only 14 percent 

when all fibers are completely disbonded. The 
significance in these results Is that the final 
value for a ^22 and a Q .33 virtually equal to 
a m . This would tend to imply that the sponge 
effect discussed earlier has little effect on 
thermal expansion. 

It is significant to note that otQj i , ag, 22 * 
and ag ,33 all approach ot^ when the composite 
is approaching a quasi -monol i thi c state. 

High Temperature (1500 °F) Ply Property With All 
Fibers Disbonding 

Consider now high temperature (1500 °F> con- 
ditions while disbonding nine fibers. E<m is 
nonlinear and degrades rapidly during initial dis- 
bonding but tends to stabilize as disbonding 
approaches 100 percent (Fig 16). This curve is 
similar in shape to that of room temperature Ejh \ 
except it is shifted down the ordinate. The total 
decrease in Egj ] is 85.6 percent and its value 
at 100 percent disbonding 1 s lower than that of 
the matrix alone. Again, this is probably due to 
the sponge effect. The transverse modulus Eg ,22 
degrades by 20.6 percent at 100 percent fiber 
disbond (Fig 16). 

Figure 17 shows the effect of fiber disbond- 
ing on the Poison's ratios at high temperature for 
all fibers disbonding. As with the room temper- 
ature results, vqj 2 increases initially and then 
begins to decrease at about 25 percent disbonding. 
The total decrease in ^ rom 0 percent dis- 

bonding to 100 percent disbonding, is 9.1 percent. 
V Q,23 decreases linearly with fiber disbonding as 
seen in Fig. 17. These material properties 
decreased by 29.6 percent when 100 percent dis- 
bonding i s reached. 

Figure 18 shows a linear decrease in Gqj 2 
with fiber disbonding. 2 decreased by 
35.7 percent when the composite is at 100 percent 
disbonding. This is the same trend seen at room 
temperature. The degradation of G& 23 , as seen 
in Fig. 18, is slightly nonlinear both at high 
temperature and at room temperature (Fig. 14). 

G £23 shows a slight increase in slope after about 
25 percent disbonding. The total decrease in 
G &23 is 34.2 percent at 100 percent disbond. 

Figure 19 depicts the change in the composite's 
thermal expansion coefficients through the range 
of 0 to 100 percent disbonding at 1500 °F. As the 
percent disbond increased, a significant increase 
in olqji is evident. At 100 percent disbonding, 
agji has increased by 799 percent and is virtually 
the same as 0^22 and ajm at 100 percent dis- 
bonding. The degradation In <xq, 22 appears to be 
minimal. The three thermal expansion coefficients 
approach the value of the matrix expansion coef- 
ficient at 100 percent fiber disbonding. 


Conclusions 


From the results of this analysis it can be 
seen that, in general, single fiber fracture and 
disbonding have little effect on most of the 
material properties. Of the material properties 
considered, the longitudinal thermal expansion 
coefficient agj 1 and the longitudinal modulus 
Eail appear to be more sensitive to fiber dis- 
bonding. This sensitivity of agj ] makes it a 
good indicator of the level of disbonding in a 
composite. If agj 1 is above predicted values, 
then the suspicion of the presence of disbonding 
of fractural fibers may be justified. The degree 
of disbonding may be approximated by comparing 
with the results shown here. 

High temperature results indicate a shift in 
the value of the material properties while main- 
taining the same rate of degradation as in the 
room temperature results. 

When all the fibers are fractured and dis- 
bonded, more significant changes in the material's 
ply properties are observed. The longitudinal 
modulus, Efjji, and vqj 2 both show nonlinear 
degradation at these higher levels of fiber dis- 
bonding. Again the sensitivity of a^n to the 
disbonding of fracture of all fibers reinforces 
the idea that it is a good indicator of the degree 
of disbonding and fiber fracture in a composite. 

This work has shown that for small amounts 
of partial bonding on fractured fibers, composite 
material properties are not significantly 
affected. However, when most of the fibers are 
fractured and disbonded, substantial degradation 
of the ply properties can be expected. 

References 

1. Caruso, J.J., "Application of Finite Element 
Substructuring to Composite Micromechanics," 
NASA TM-83729, 1984. 

2. Caruso, J.J., and Chamis, C.C., "Superelement 
Methods Applications to Micromechanics of 
High Temperature Metal Matrix Composites," 
AIAA 29th Structures, Structural Dynamics, 
and Mater i a 1 s , Conference, Part 3, AIAA, 

1988, pp. 1388-1400. 

3. Gockel , M.A., ed., MSC/NASTRAN Handbook for 

Superelement Analysis: MSC/NASTRAN Version 

6J_, The MacNeal -Schwendl er Corporation, 

Los Angeles, 1982, p. 1.1-1. 

4. Hopkins, D.A., and Chamis, C.C., "A Unique 
Set of Micromechanics Equations for High 
Temperature Metal Matrix Composites," NASA 
TM-87154, 1985. 


I 


5 


TABLE I. - COMPOSITE MATERIAL 
PROPERTIES NO DISBONDING - 
ALL FIBERS LOADED 
P100-GRAPHITE/ 

COPPER 


TABLE II. - ROOM AND HIGH TEMPERATURE CONSTITUTIVE MATERIAL 
PROPERTIES PI OO-GRAPHITE /COPPER 
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Figure 3. Center fiber disbonding showing 2.78% disbonding 
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MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
MODULUS - Elll, E122, E133 



PERCENT DISBOND 

Figure 4. Effects of center fiber disbonding at room temperature 
on moduli 
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POISSONS RATIO 


EFFECT OF FIBER DEBONDING ON 
POISSONS RATIO - i412, i413, 1/123, i432 
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Figure 5. Effects of center fiber disbonding at room temperature 
on Poisson's ratio 





SHEAR MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
SHEAR MODULUS - G121, G131, G123 



PERCENT DISBOND 

Figure 6. Effects of center fiber disbonding at room temperature 
on shear moduli 
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TEC (xIO - 8 in/in/°F) 


EFFECT OF FIBER DEBONDING ON 
THERMAL EXPANSION COEFFICIENT (TEC) - alll, cd22, al33 



PERCENT DISBOND 


Figure 7. Effects of center fiber disbonding at room temperature 
on thermal expansion coefficients 
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EFFECT OF FIBER DEBONDING ON 
MODULUS - Elll, E122, E133 
HIGH TEMPERATURE (1500 °F) 
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Figure 8. 


Effects of center fiber disbonding at high temperature 
on moduli 
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POISSONS RATIO 


EFFECT OF FIBER DEBONDING ON 
POISSONS RATIO - i412, i/113, i/123, i/132 
HIGH TEMPERATURE (1500 °F) 



PERCENT DISBOND 


o 

(D 

o 

o 


e 

o 

O 

o 



1 / 112 . 

]A13 

1/123 

i/132 


Figure 9. 


Effects of center fiber disbonding at high temperature 
on Poisson's ratio 
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SHEAR MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
SHEAR MODULUS - G121, G131, G123 
HIGH TEMPERATURE (1500 °F) 



PERCENT DISBOND 


Figure 10. Effects of center fiber disbonding at high temperature 
on shear moduli 
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TEC (xIO -6 ln/In/°F) 
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Figure 11. Effects of cepter fiber disbonding at high temperature 
on thermal expansion coefficients 
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Figure 12. Effects of nine fiber disbonding at room temperature 
of moduli 
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Figure 13 . Effects of nine fiber disbonding at room temperature 
on Poisson’s ratio 
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Figure 14. Effects of nine fiber disbonding at room temperature 
on shear moduli 
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Figure 15. Effects of nine fiber disbonding at room temperature 
on thermal expansion coefficients 
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Figure 16. Effects of nine fiber disbonding at high temperature 
on moduli 
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Figure 17. Effects of nine fiber disbonding at high temperature 
on Poisson's ratio 
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Figure 18. Effects of nine fiber disbonding at high temperature 
on shear moduli 


22 




(do/ u !/ u ! 9-OI-X) oa± 


EFFECT OF FIBER DEBONDING ON 
THERMAL EXPANSION COEFFICIENT (TEC) - cdll, otl22, «133 
HIGH TEMPERATURE (1500 °F) 



PERCENT DISBOND 

Figure 19. Effects of nine fiber disbonding at high temperature 
on thermal expansion coefficients 
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